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We discuss conformal field theories (CFTs) in rectangular geometries, and develop a formalism 
that involves a conformal boundary state for the 1 + Id open system. We focus on the case of 
homogeneous boundary conditions (no insertion of a boundary condition changing operator), for 
which we derive an explicit expression of the associated boundary state, valid for any arbitrary 
CFT. We check the validity of our solution, comparing it with known results for partition functions, 
numerical simulations of lattice discretizations, and coherent state expressions for free theories. 



I. INTRODUCTION 

Boundary conformal field theory (BCFT) is a subject 
whose importance has grown over the years, both on the 
formal and on the applied side. It does for instance play a 
fundamental role in the axiomatization of conformal field 
theory (CFT) [1], in our growing understanding of loga- 
rithmic conformal field theory [2, 3], or in the relation- 
ship between conformal field theory and the Schramm 
Loewner Evolution formalism [4]. It is also central to 
our understanding of the Kondo effect [5] , of the physics 
of quantum impurities or the Fermi edge singularity [6] , 
and, more recently, of local and global quenches in one 
dimensional quantum systems [7, 8]. 

Two geometries are most naturally used in BCFT. One 
of them — more natural from the point of view of Euclid- 
ian field theory or statistical mechanics — is simply the 
upper half plane, with the theory defined e. g. only for 
Im z > 0, together with some boundary conditions at 
Im z ~ 0. The other involves a cylinder (sometimes con- 
sidered instead as an annulus), and typically describes 
the physics of an open ID quantum system at finite tem- 
perature. 

Of course, there are variants of these geometries. For 
instance, the physics of the cylinder can be described 
with an (imaginary) time evolution along the axis. In this 
case, the ID quantum system is closed, and the presence 
of the boundary is encoded in a boundary state. Such 
states have given rise to many developments, in the con- 
text of efforts to classify all conformal boundary condi- 
tions in particular. 

We shall be concerned in this paper mostly with sys- 
tems on a rectangle, which thus have trivial topology, 
sharp corners, and, in general, four different boundary 
conditions on their four edges. While in principle this 
situation can be tackled via conformal mappings of the 
half plane (see below), few of its features have actually 
been studied in detail. The situation is, on the other 
hand, clearly interesting in its own right. The rectan- 
gle geometry is, for instance, the natural one to consider 
in the case of quenches for II? quantum systems with 



open boundaries. In the 2D point of view, it is the sim- 
plest geometry to study transport properties of network 
models for Anderson localization. More fundamentally, 
the rectangle provides a natural way to study and in- 
terpret the conformal blocks of four point functions by 
inserting four different fields at the four corners. This 
is particularly useful for instance to connect geometrical 
correlators in the Self Avoiding Walk (SAW) problem to 
CFT, in particular in the logarithmic (indecomposable) 
case. 

We will focus on these and other geometrical features 
in subsequent and rather technical work to appear soon. 
The present paper is devoted to the exploration of the 
simplest — and yet very rich — aspects of this problem re- 
lated with the boundary state for a theory defined on a 
segment. 

II. GLUING CONDITION 

For a CFT in the complex plane C, one can implement 
a boundary with the requirement that there is no energy 
flow across it. In other words, one of the components of 
the stress-tensor vanishes along the boundary: Tj|j_ = 0. 
In holomorphic/anti-holomorphic components this gives 
the constraint T(z) = T(z) on the boundary. In radial 
quantization, a circular boundary at radius \z\ = 1 can 
be encoded in the form of a boundary state \B P ). Here 
the subscript p stands for periodic since the boundary is 
defined on a circle. The constraint is 

(L n -L_ n )\BP) = 0, (1) 

which is usually referred to as the conformal invariance 
of the boundary condition, or gluing condition, as it glues 
the modes of the chiral part of the CFT with the anti- 
chiral ones. A basis of solutions of the linear system of 
equations (1) is given by the so-called Ishibashi states, 
which arc particular combinations of left-right symmet- 
ric Virasoro descendants of primary fields. Of course, 
to obtain the allowed boundary states \B P ) themselves, 
more conditions have to be implemented (see for example 
[1, 5] for reviews). 
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In this paper, we shall extend the formalism of bound- 
ary states to the case of an open system. Namely, we 
are now considering the CFT on a strip of width L, and 
we want to find a state in the theory that encodes the 
boundary perpendicular to the direction of the strip (see 
Fig. 1). In the 1 + 1 Hamiltonian description of the CFT, 
the boundary state now lives in the Hilbert space of the 
Id theory defined on a segment instead of a circle. To pro- 
ceed, we choose coordinates z = x + iy in the plane and 
consider the semi-rectangular region < x < L,y > 
of Fig. 1, and derive the gluing condition for the stress- 
energy tensor in this geometry (for analogous discussions 
see [9-11]). Conformal invariance of the boundary im- 




FIG. 1. The semi-rectangular geometry denning the boundary 
state. 

plies then T xy = on every side. This is realized at the 
boundaries z = iy and z = L + iy, y > 0, by setting 



T{z) = T(2L - z) . 



(2) 



and imposing periodic boundary conditions in the x di- 
rection of period 2L; T(z + 2L) = T(z) (same for f). 

Care is needed to treat properly the effect of corners. 
The stress tensor T has a singularity as its argument ap- 
proaches the corners. To see this we start more generally 
by the singularity in the upper half plane when an oper- 
ator of weight h is inserted at the origin. In this case the 
most singular term is 



T(w) 



(3) 



We then fold the upper half plane by the mapping z = 
w 1 / 2 to have a corner in z — 0. After using the transfor- 
mation law of the stress tensor T(z) — T(w)(dw/dz) 2 + 
(c/12){u>; z}, where {w; z} is the Schwarzian derivative 



{w;z} 



d 3 w/dz 3 
dw/dz 



d 2 w/dz' 
dw/dz 



(4) 



As a consequence the condition on T defining the 
boundary state \B°) at y = is 



; T{x) - T(~x) + Airi ( hiS'(x) + h r S'(x - L) ) ) \B°) = , 

(6) 



where we have used 
1 1 



Aiex 



(x — it) 2 (x + ie) 2 (x 2 + e 2 )' 



-> -2mS'{x) , (7) 



and defined the "effective conformal weight" at the cor 
ners 

c 



h:=2h- 



16 ' 



(8) 



hi (h r ) being the weight of the operator inserted at the 
left (right) corner. Since the system is periodic with pe- 
riod 2L, we can go to a mode expansion 



T(x) 



J2 L n er mx/L 



to get the gluing condition (for n £ Z >0 ) 



L_ n -2n{h l + (-) n h r ))\B°) = 



(9) 



(10) 



Note in particular that in this case one identifies not 
only L n with L_„ but also L n with L_„, since the semi- 
rectangular geometry is obtained by two consecutive fold- 
ing of the plane. 

Calling the boundary conditions on the sides of the 
semi-rectangle a, &, c, and the boundary condition chang- 
ing operators sitting at the corners 4>f{0) and tf W {L) 
(with opportune labels i,j), the boundary state \B°) 
will gcnerically live in a direct sum of vector spaces de- 
termined by the fusion tfi^ b (g) </>^' c . In this paper we 
consider the case of homogeneous boundary conditions 
a = b = c, where the identity operator sits at the corners 
(hi = h r = 0), so that equation (10) becomes 



- [! + (-!)"] HB°>=0, 



(11) 



we get 



to be solved within the Verma module of the vacuum |0) 
of the CFT. The solution of the homogeneous case will 
then be of the form G |0), with G a certain expression in 
terms of the Virasoro generators L_ n 's. When boundary 
condition changing operators are present, the boundary 
state will simply be G <j)^ b \w\)^ c (11)2) |0), where the 
two points w\ and W2 he on the boundary and correspond 
to the images of the two corners at z = and z = L under 
a conformal mapping to the upper half plane. Further 
details about this general situation will be discussed in a 
sequel [14]. 



T(z) 



Ah 



(5) 



Also if h — 0, at a corner there is an anomaly, which re- 
flects itself in a non-trivial scaling dependence of physical 
quantities [12, 13]. 



III. BOUNDARY STATE IN BOUNDARY CFT, 
AND THE RECTANGULAR BOTTOM 



In this section we give the explicit form of the bound- 
ary state I B°) which solves the constraint (11), following 



3 



GW I pi. _ 




FIG. 2. A deformation of the boundary around the origin can 
be encoded in a conformal mapping g : H' — >• H. 



the discussion in [15]. 
instance, in [4, 16]. 



Similar tricks have been used, for 



A. CFT in the half-plane, deformations of the 
boundary around the origin 

Let us consider some boundary CFT defined in the 
half-plane H. We make use of radial quantization in H. 
The vacuum of the theory can be written formally as a 
path integral 



|0> 



[#(1*1 < 1,2! G 



I0(M = 1)>, (12) 



where 4>{\z\ = 1) stands for the configurations of the fields 
of the theory on a semi-circle centered at the origin and 
of radius 1 . The weights of the different configurations is 
given by a Gibbs distribution e~ s ^'. In this section, we 
consider deformations of the half-plane around the origin, 
such as the one shown in Fig. 2. More precisely, we take 
some domain H' = H \ K with K a domain included in 
the semi-disc of radius 1 centered on the origin, such that 
H' is simply connected. Our goal is to find an expression 
for the new state 



') = / [d<j>{\z\<l,z£ 



\m = i)) • as) 



In this paper we do not consider insertions of boundary 
condition changing operators or other operators in the 
bulk. Therefore, we expect |H') to be given by some 
linear combination of the vacuum |0) and its descendants 
only. In other words, there must exist some operator Gw 
built out of the Virasoro generators such that 



Gw |0> 



(14) 



Our goal is to find how to construct Gw ■ 

By Riemann's mapping theorem, there is a conformal 
mapping from H' onto H. If one requires that the behav- 
ior of g(z) as z — > oo be 



Ol u>3 

g(z) =z + — + — + — + 



02 
2 



(15) 



then the mapping g is unique. 

Now, let us assume that we have a continuous family 
of deformations of the half-plane {H t }o<t<i, such that 



Ho = H and Hi = H'. For each t there is a mapping g t 
from H t onto H with the above asymptotic behavior. The 
composition g^ 1 o g t+dt gives us an infinitesimal mapping 
from Ht+dt onto H t , which can be expanded as 



1 ° 9t+dt{z) = z + + + 



(16) 



This allows us to relate |Ht) and |Ht+dt), using the defi- 
nition of the stress-tensor in boundary CFT, namely the 
variation of the action S[4>] under a small transformation 
z i ^ z ~\~ a(z) (see e.g. [17]) 



SS = ^-<j>a{z)T{z)dz- 



(17) 



which gives for a(z) — dt (b\/z + b 2 /z 2 . . .) 



\m +dt ) = / [#(|zi < l,z G H t )](l - SS)e~ s |0(|*| = 1)) 
= |H t ) 

\J> a(z)T(z)dz - j> a(z)T(z)dzj |H t ) 

= \W t ) - dt (&!.L_a + & 2 i- 3 + & 3 i-4 + • ■ • ) |H t ) . 
We arrive at the differential equation 

^ \Mt) = - (6iL_ 2 + 6 2 L_ 3 + b 3 L_ 4 + ...) \W t ) (18) 
dt 

or, in terms of the operator Gt = Gn t 

(&!L_ 2 + 6 2 i_ 3 + &3^-4 + ...)G t . (19) 



dt 



Note that (16) and (19) provide us with a differential 
equation which, at least formally, allows us to find the 
operator Gw = Gt=i- This differential equation will play 
an important role in the next sections. 



B. The boundary state 

In general, it is not possible to find an explicit solution 
to the above differential equation given some geometry 
H'. For certain cases, however, solving the equation is 
easy. 

Let us start with the case of a vertical slit of height \p2± 
in the half-plane, namely the domain H \ [0, iy/2t\ , see 
Fig. 3. It can be mapped onto the half-plane by gt : z H> 
Vz 2 + 2t. Note that gt{z) = z + - + . . . as z — > oo, so 
it has the required asymptotic behavior. Differentiating 
this function with respect to t, we get 



9i 



1 ° gt+dt{z) = z + 



dt 



(20) 



According to the previous section, this corresponds to an 
operator Gt which is given by 



dt 



Gt 



-L-vGt 



(21) 



FIG. 3. The half-plane minus a slit of height i/2t, namely 
H \ [0, i^/2t\ , is the simplest example of a domain which can 
be encoded in a boundary state: it corresponds to e~ tL ~ 2 |0). 

with the initial condition Gt=o = 1. Wc thus find Gt — 
e~ tL ~ 2 . In particular, the state corresponding to a slit of 
height 1 is e~? L - 2 . 

It turns out that one can easily extend this trick to a 
larger set of domains with slits. Let us define the half- 
plane minus k — 1 slits of size 2 1 / fc 

e fc = e\{z|z*e [-2,2]} . (22) 

See Fig. 4 for an illustration of the multi-slit geometry, 
and note that Hi = H. The point of this definition is the 
following. We define the functions 

g k {z)= (z k + 2) 1/k (23) 

and wc make the observation that g 2 N is a confor- 
mal mapping from H 2 n onto H 2 w-i. Moreover, these 
mappings have the correct asymptotic behavior at in- 
finity, and one can build a family of mappings g k ,t = 

(z fc + 2t) 1 ^ k which gives rise to 

9ll ° 9k,t+dt (z)=z+- p-y (24) 

and therefore to a differential equation for the operator 

Gk,t = Gu k _ t 

-^G M = ~L_ k G k>t (25) 
at k 

which gives in the end e~^ L ~ k . 

As already mentioned, g 2 N is a mapping from H 2 n 
onto H 2 n-i, so by composing several of those, we see 
that g 2 54 • ■ ■ <?2 N is a mapping from H 2 n onto the 
half-plane H. Thus, it is now straightforward to write 
down the boundary state |H 2 iv) with 2^ — 1 slits of size 

2 l/2« 

|H 2 «) = e -^hr L -2" . . . e -^£-4 e -L-2 | ) _ (2 6 ) 

The connection between these states and the main 
topic of this paper appears when one considers the limit 



FIG. 4. The domain Hg is the half-plane minus 7 = 8 — 1 slits 
of size 2 1/s . 




-11 -2 2 



FIG. 5. Mapping of semicircular region T> to the upper half 
plane H f{z) = z + z' 1 



N — >• 00. Indeed, one can check that 
g 2 a 54 o • • • o 32" (z) = 




when \z\ — > 00. We therefore have lim 32 <?4 • • • 

A^— ¥00 

92 N ( Z ) = z + z^ 1 when \z\ > 1. The function z i-> z + z _1 
is a conformal mapping from V := H \ {«, |z|<l} to 
H, see Fig. 5. It is thus tempting to conclude that the 
boundary state \B°) which corresponds to this domain 
(Fig. 6) is [15] 

\B°) = lim e -w=T L -2" . . , e -^ L - i e^ L - 2 10) . (28) 

N—*oo 

At first sight, this expression might look ill-defined, how- 
ever it is not. Indeed, it is obvious that, expanding the 
different exponentials in powers of the Virasoro modes, 
one gets a finite sum at each level. For instance, our 
boundary state is, up to level 4 

\B°) = |0) - L_ 2 |0) - h^\0) + ii 2 _ 2 |0) + . . . . (29) 

The foregoing argument leading to cq. (28) is of course 
not entirely satisfactory, since it relies on an a priori un- 
comfortable identification of the half disk with an infinity 
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FIG. 6. The boundary state \B°) is obtained in the N 
limit. 



of slits, supported by the asymptotic behavior (27). This 
can hold only in "some sense" — we believe, in the sense 
that all correlation functions in the CFT for the geome- 
try of interest can be obtained using eq. (28). We explore 
this further in what follows. 



C. Deformed gluing conditions 

In this section we will derive the gluing conditions sat- 
isfied by the stress tensor in the multi-slit geometries H 2 w 
introduced in section IIIB. This will lead to a proof of 
the fact that the boundary state (28) satisfies the gluing 
condition for the semi- rectangular geometry, eq. (11). 

The region H 2 « can be mapped to the upper half plane 
by the conformal mapping (27), whose expansion at z — > 
co is 



w = Jn(z) ■= z H h V 

z 4 z 

p=K 



Ok 
k 



(30) 



where we have introduced K := 2 N+1 — 1. Its inverse has 
the following expansion at w — > oo: 



h N (w) 



w 



+ Vi 



y §h. 

p—K 



(31) 



For the general argument given below we will not need 
the explicit expression of 's and /3 k 's so we leave them 
unspecified. We note however that closed expressions of 
the above mappings can be found: 



f N (z) = 2 cos 2 



,-JV 



arccos 



(32) 

h N (w) = (2 cos (2 N arccos (f ))) 2 " , (33) 

which eventually allow to write explicitly the coefficients 
appearing above and in the next computations. Let us 
now derive the gluing condition for the multi-slit geom- 
etry by using a slightly different point of view from the 
one of section II. 

If we call X = ■ ■ ■ 4>{ z ii z\) ■ ■ ■ tp{x) ■ ■ ■ a chain of ar- 
bitrary bulk and boundary operators, by definition our 
boundary state (26) is such that 



x|]&v) = (o|l|o) , 



(34) 



where X = G m n XGm 2N is the conjugation by the opera- 
tors implementing the conformal mapping. Note that the 



out- vacuum (0| is invariant under this mapping because 
In(z) ~ z when \z\ — > oo. Inserting the stress tensor in 
this correlator we have 

(0\XT(z)\U 2 n) 

= (ti N (w)r 2 [(0| XT{w) |0) - ^{h N ; w} (0| X |H 2 «) 

(35) 

where {/iat;w;} is the Schwarzian derivative (4). We use 
this starting point to compute 



<0| X (L n - L_ n ) 



= (01 X 



X 



dz 

— z(z n - z- n )T(z) |H 2 ») 
2ni 

dw 



2iri 



J N . n {w)T(w)\0) 



c f dw 

12 t- ^~ Jn ^ w ^9n; w} (0| X 



(36) 
(37) 

(38) 

(39) 



The contour C above encircles the points w = ±2 and lies 
in the region of validity of the Laurent expansion (31), 
and we have defined 



J 



N.. 



N 



'N 



Further one can check that 

J N ,nH = (™ 2 ~ 4)^-1 (I 
where 

i 

/ n — 1 — k 



wP 

p—K — n 



(40) 



(41) 



t4-i f 



k 



i-l-2fc 



(42) 



is the Chebyshev polynomial of the second kind and a' p 
are opportune coefficients. Note then that as N — >• oo 
the function in eq. (41) reduces to a polynomial of degree 
n+ 1. We focus now on the the integral in equation (38). 
Computing the residue at the origin, the polynomial part 
of Jn,u gives zero, while mapping the contribution of 
the other non-analytic part back to the original multi- 
slit geometry one finds 

{0\X f J^J N , n (w)T(w)\Q) 

OO 

= A*<0| *£- P -i |H 3 *) (43) 

p—K — n 



C 

12 



Y, a ' P f ^- p {^;^}(0|X|H 2 ») , 

/ s~i Z7TZ 



p—K—n 

where // p 's are some coefficients whose explicit expression 
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is not needed. Plugging this back in (38) we have 

oo 

(0| X (L n - L_ n ) \W 2 n) = M0|^£- P -i|H 2 

p— K —n 

iW ,(w 2 ~ 4)C/„_ 1 (f ) {h N ; w} (0| X |H 2 ») 



12 Jc 2iri 



(44) 



We now separate the N — > oo result from the finite N 
contributions in the Schwarzian derivative to get an ex- 
pression of the following form: 



{9n;w} 



6 



-+ V 



(45) 



Both terms in the equation above give a non-vanishing 
contribution: 



(w 2 ~ 4)C/„_ X (f ) {h N ; w} (0| X |H 2 ») 



dw 
! Jc 2ttI 



(1 + (-!)«) + A„, 



(46) 



where we called A n the number resulting from the the 
integration of the second term in (45). Clearly one has 
A n = if n < K. In the end one finds the following 
gluing condition for the multi-slit geometry: 



(i n -L_ n + y [l + (-l) n ]) |H 2 «) 

OO \ 

^ /ipi_p_i + A n \M 2 n) . 
ip—K — n J 



(47) 



This deformed constraint reduces to the one of eq. (11) 
for the semi-rectangle when N — > oo. This is enough to 
prove that the boundary state (28) obtained as limit of 
the geometries with 2 N — 1 slits indeed satisfies the gluing 
condition (11). 



D. Comparison with partition functions 

A next obvious step is to consider the amplitude as- 
sociated to the boundary state eq. (28) defined in the 
previous section and compare it with the known result 
for the universal part of the partition function on a rect- 
angle. We recall first this result. Since we consider the 
case of the same boundary condition on each side, we 
note that the partition function for a rectangle of length 
L and width L' will have to be a modular form. If we call 
the universal part of the partition function Z-jz{L, L') 



e hLL' e f*(L+L')Z K (L,L'), 



(48) 



where /b and / s are bulk and surface energies, we have 
Z n (L,L') = Z n (L',L). (49) 



Define as usual r := iL'/L, q 



; 27 ™ T and the 



Dedekind eta function as r?(r) = q 1 / 24 l\™ =1 (1- q n ). The 
partition function Zn(L, L') was computed in [18] and up 
to possible proportionality coefficients is: 

-c/2 



Z n {L,L') = L c 'S{T) 



(50) 



c(c + 6) 2 
8 q 



As it was observed in [19], this partition function can be 
also derived directly from modularity arguments. Indeed, 
on top of invariance under modular inversion, as we have 
discussed the presence of corners corresponds to an effec- 
tive weight c/16 for each corner, which fix the non-trivial 
scaling of the partition function. 

Now we can take the scalar product of the boundary 
state with itself to form the amplitude 



B \q L °- c ' 24 \B°) (51) 
c n q n , (52) 



where q := yjq is the relevant combination of L' /L ap- 
pearing in the transfer matrix on a strip, and the arrow 
is the direction of imaginary time. The following rela- 
tion between this amplitude and the partition function 
Z n (L,L') should hold 



A n (r) = L-^Zn^L') = V (ry c / 2 . 



(53) 



We have verified this relation by computing the first co- 
efficients c n in (52) using the commutation relations of 
the Virasoro algebra and comparing them with the power 
series of r]~ c ^ 2 , confirming the validity of our derivation 
up to a very high order (n — 52) in q. We note that as 
expected from the left-right symmetry of the boundary 
conditions of the problem, the boundary states couple 
only to descendants of the identity of even level (so that 
c 2 n+i = in eq. (52)). 



E. Convergence and multi-slit geometry 

In this subsection we discuss further the N — > oo limit 
illustrated in Fig. 6, by computing amplitudes involving 
the Unitized boundary state 



I 2 «) = e'^t ~* N - 



corresponding to 2 N — 1 slits. We are interested in the 
generating functions 

P N (q) = ((M 2 «|q L «/ 2 |H 2 «)) = $>£V ■ (55) 

k=0 

Since the expected limit is 

X . oo 

^lim P N (q) = q^Mr)- 1 = R — - = J> 



m— 1 



k=0 



(56) 
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where pk is the number of partitions of the integer k, we 
expect that the p^ will somehow converge to pk- We 
shall see below how this occurs. 

The case N = 1 is easily dealt with analytically. Using 
induction one can prove that 



k-i 



\L k 2 L k _ 2 \0) = ¥ l[(8p + c) 



(57) 



p=0 



Developing the exponential e L ~ 2 one then shows that 



P 1 (q) = (l-4q)- 1 ^ = l + q+-q 2 + ... . 

In the next case, N = 2, a direct computation gives 

33 

P 2 {q) = 1 + q + 2q 2 + 3q 3 + —q 4 + ... . 



(58) 



(59) 



By automatizing the computations in the Virasoro alge- 
bra using Mathematica we have been able to obtain the 
P/v(<?) up to order q 26 . In the case of £2(9) the results 
are consistent with the conjecture 



(l + 2g) 1 / 2 (l + 4 g 2 ) 5 / 8 
(1 - 16g 4 ) 3 / 4 



For higher N we have not been able to conjecture — let 
alone derive — such exact expressions. We find 



P 3 (q) = 1 + q + 2q 2 + 3g 3 + 5g 4 + 7q 5 + llq 6 



7 245 s 
15<Z 7 + — q 8 



14988511 



and 



Pkq 



k=Q 



, 4005 16 
16 q 



512 



27657 



-q 



.20 



26 



.(61) 



.., (62) 



while P${q) agrees with Y^h=oPkq a ^ l eas t up to order 
q 26 . 

So we observe from these examples that the coefficients 
p~^ are non-negative rationals that coincide with the 
integers pk for k = 0, 1, . . . , 2 N — 1. The first deviating 
coefficient occurs for ko — 2 N , and we have then p^' > 
Pk - We conjecture that these observed properties hold 
true for any N > 1. 



IV. THE CASE OF FREE THEORIES 

While the expression (28) is general, there are more 
natural ways to think of the boundary state in the case 
of free theories. There, like when the boundary is a circle 
instead of a segment, expressions as coherent states over 
the bosonic/fermionic modes are possible. The compari- 
son with results in the first section is intricate, and leads 
to remarkable identities. We note that some of the co- 
herent states expressions presented in this section have 
an important overlap with [9-11] [20]. 



A. Free boson 

The simplest case one can deal with is the free boson. 
We introduce then harmonic oscillators a„ , n € Z sat- 
isfying [a m ,a n ] = mS m+ni0 , a m \0) = if m > 0, and 
ao|0) = 0, with which we represent as usual the Virasoro 
generators: 

L n = — y ^ a n — m a m , if n 7^ , 
La = — a 2 , + ' S y^ j a,- m a m . 



(63) 



Substituting this expression in (10), we realize that the 
solution should be of the form: 



|B;)=exp(-^i- a 2 _„ 

\ n>0 / 



(64) 



To verify this one simply has to use twice the following 
relation 



{a m + o_ m ) \B%) = , m>0. 



(65) 

(60) The extra term n/8(l + (-)") comes in because if n = 2p 



we have a 2 \B°>) 



p)\BT)- Another straightfor- 



ward exercise is computing the amplitude associated with 
this boundary state. We have 



{B o\qL -l/2A wo) = 



1/48 



n 

m>0 



(2s)! (q 



(66) 



in agreement with (53). Alternatively one could also ex- 
press L n in (28) using (63) and check that the two ex- 
pressions agrees for c = 1 . 

Further we note that this result for the partition func- 
tion of the free boson can also be obtained by computing 
the determinant of the Laplacian with, say, Neumann 
boundary conditions on all four sides. To our knowledge 
this result appeared first in [21] (there is factor of LL' 
there due to the way the zero mode is subtracted): 



det(-A) = L- 1/2 77(r) 



(67) 



We note the explicit presence of the anomaly term 
L c / 4 in the partition function given by Zn(L,L') = 
(dct(-A))- 1 / 2 . 

It is interesting to observe that eq. (67) corresponds, 
using the matrix-tree theorem, to the partition function 
of spanning trees, equivalent in turn to dense polymers or 
symplectic fermions — all variants of c = —2 CFT. (More 
generally, partition functions for c = — 2 with different 
boundary conditions on the sides will be given by taking 
the power —2 of the expressions for the free boson, since 
the boundary condition changing operators have dimen- 
sions 1/16 for c — 1 and —1/8 for c = —2. This will be 
exploited in a further paper [14].) 
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B. Majorana fermions 

We consider now the case of Majorana fermions. We 
will go carefully through the derivation of the gluing con- 
dition which turns out to be more tricky in this case, 
and then derive the boundary state by computing the 
correlator in the rectangular geometry using a conformal 
mapping to the upper half plane. 

First, we fix notations and conventions. The fermionic 
modes are defined via the expansion of the bulk fermionic 
fields 

rel r£l 

where I = Z in the Ramond (R) sector and I = Z + 1/2 
in the Neveu- Schwartz (NS) sector. The OPEs 

ip(z)ip(w) = , ijj(z)'ijj(w) = (69) 

z — w z — w 

lead to 

{lp r ,1p a } = {Vv,^ s } = S r+St0 - (70) 

Let us recall first the condition in the cylinder geometry 
with Dirichlet (D) or Neumann (N) boundary conditions 
on the bottom [22]. For simplicity we restrict to the 
NS case, where fermions are periodic in the plane, and 
thus antiperiodic on the cylinder. The cylinder can be 
obtained from the plane (with complex coordinate z) by 
the mapping w = iln z or z = e~ lw . Set w = x + iy. The 
mode expansion now reads 

^(x,y) = (-i)^ 2 J2Ae irx - rv , (71) 

r 

i,(x,y)=i 1 / 2 Y,*Pre- trx - ry . (72) 

r 

The boundary conditions are then 

i>(x,y = 0) = eip(x,y = 0), (73) 

where e = 1 for N and —1 for D. The boundary state 
associated to these equations is well-known and has the 
form of (fermionic) coherent states: 

\B^) oc exp I ei 5>- P -i/ 2 V~- P -i/2 10), (74) 
V p=o J 

where 

VV+i/ 2 |O) = 0, ^+i/2|O)=0, P eZ> . (75) 

Now imagine doing the same problem in a geometry 
rotated by ±7r/2, obtained by the mapping w' = ±iw. 
In this rotation, we use the general formulas for trans- 
forming Majorana fermions, which are objects of dimen- 
sion 1/2, and thus ip(w')(dw') 1 ^ 2 = ip{w)(dw) 1 / 2 , and 
similarly for i/j. Hence the equations characterizing the 



boundary state now would read ip = ieiip instead of 
(73): when dealing with fermions, the D or N boundary 
condition is not represented by a unique equation, but 
depends on the orientation of the boundary. 

Let us go back now to the geometry of interest of Fig. 1 , 
and for definiteness choose, say, D boundary conditions 
on all sides and L = tt. The case with N boundary con- 
ditions can be treated in a similar way, and lead in the 
end to the same result. According to our discussion we 
have 

ip(x = 0, y) = i%jj(x = 0, y) , y e [-co, oo] 

xjj(x = n,y) = -it/j(x = ir,y) , yg [-00,00] (76) 

Fully open boundary conditions require identification 
(up to monodromy conditions) of left and right modes 
from both the direct and crossed channels points of view. 
It is useful therefore to define 

i/>(x,y) = -i$(-x,y) , x€[— tt,0] (77) 

so now we can re-express everything in terms of only one 
type of fermions — we choose tp. Moreover, by construc- 
tion, tp is now regular at the origin. By symmetry we 
can also extend fermions on the other side of the right 
boundary by defining 

ip(2Tr-x,y)=i'tjj(x,y) (78) 

so t/j is now regular at x = tt. Using (77) gives finally 

%j){x + 2-k) = -ip(x) , xe[-ir, tt] (79) 

so that we now have a single species of fermions defined 
on the circle [— tt,tt], with NS (antiperiodic) boundary 
conditions. 

All the manipulations so far are there to handle the 
boundary conditions at x = and x = tt. We have now, 
on top of this, to handle the boundary condition at the 
bottom of the system, that is for y = 0. Taking e = — 1 
in (73), from (77) we have 

[i>(x, y = 0)+ ii>(-x, y = 0)] \B°) - , x G [0, tt] 

(80) 

Going through the definitions carefully shows that the 
relative sign of the two fermion terms in this equation 
switches in the interval x £ [— 7r,0]. This is in line with 
the fact that fermion correlators change sign under re- 
flection. So the final equation has to be for x G [— 7T, tt] 

[iP(x, y = 0) + sign(x)^(-x, y = 0)] \B° 4 ) = . (81) 

The sign function introduces serious complications with 
respect to the case of ordinary boundary states. No- 
tice that it disappears when considering conditions for 
the Virasoro generators, since the stress-energy tensor is 
quadratic in the fermions. We represent the sign function 
through the Fourier series 

00 

sign(z)= a ^ mX (82) 
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with 



02m = 



2i 



(83) 



(2m + l)ir 

We can introduce the matrix A defined through 

I _ ryyn+n+1 
An+l/2.n+l/2 = / 1 ~ 7V~ , m + 71 + 1 ^ 

' ' ' 7r(m + n + lj 



At»+1/2,-it»-X/2 — 



(84) 



so the equation satisfied by the boundary state reads, in 
terms of the fermionic modes 

NWl/2 + X! j4 ™+l/2,«+l/2V'™+l/2^ \%) = ■ (85) 

Now decompose the A matrix (whose labels are half odd 
integers) into four blocks corresponding to the label signs: 



.4 



&m,n Urn. Tl 



(86) 



where we have defined blocks a m n = ^4 m +i/2,n+i/2i an d 
6 TOj „ = ^4 m +i/2,-, l -i/2i with to, n > 0. Take now to > 
(for to < one gets an equivalent condition) and consider 
eq. (85), which we can rewrite as 

V-Wl/2 + ^ a m ni>n+l/2 + Kinlp-n-1/2 \B%) = . 
n>0 J 

Collecting the positive and negative fermion modes and 
inverting the system gives the equation 

vwi/2 + E [( 1+a r l6 L„^-«-i/2 |£ ^> = - 
n>o y 

Now set G := —(1 + a)~ 1 b = —G T . The solution of this 
system is (up to proportionality coefficients) the coherent 
state 

\B$ = exp [ Gm,n^-m-l/2^-n-l/2 I |0) . (87) 

\0<m<n J 

With this approach the matrix G m ^ n is somewhat un- 
determined. We present then a different strategy to com- 
pute it using correlation functions. As in section III we 
define the boundary state in the geometry denoted by T> 
(fig. 5 (left)), and consider the mapping to the upper half 
plane H (fig. 5 (right)), viz. w = z + z" 1 . We thus con- 
sider D boundary conditions on all sides of the boundary 
of T>, so that \B3,) should be the state representing this 
boundary. Calling \0) the ground state of the model at 
infinity in H we have 



/ dwi dw 2 \ X ^ 2 
\ dzi dz 2 



(88) 

(i}(w 1 )iP(w 2 ))a , (89) 



where on the rhs the correlator is evaluated in H with the 
usual D boundary conditions. Using that this correlator 
must simply be l/(w'i — w 2 ) (since it involves only right 
movers, which are not affected by the boundary), and 
evaluating derivatives gives straightforwardly that 



(OM*iM**)|B$) 



1 



1-^/1- 



1 

21 22 



Z1—Z2 1 

From the result (87) we can write the boundary state as 
dz f dz' 



exp 



2in T 2iir 



-i>(z)G(z,z'Mz') : \Q) 



_ (90) 

where G(z, z') corresponds to the generating function of 
the numbers G m ^ n introduced above. The left hand side 
of cq. (88) can now be evaluated using Wick's theorem. 
We find in the end that 



G(zi,z 2 ) 



1 



'1 - 



1 - 



2{z 2 - zi) 



f 



1 



- 1 



(91) 



We need to evaluate this expression in the domain 
|zi|, |z 2 | > 1 (because of the geometry of our problem), 
and expand it as 



G(zi 7 z 2 ) 



00 „ 
2z x z 2 ^ zV l z™ ' 

1 z m.n=0 1 2 



(92) 



(Observe that on top of G mn = —G nm , we have G mn = 
if m + n is even.) As a result, we now have the generating 
function for the the quadratic form appearing in (87). 
The first few values of G m ^ n read: 



G 



03 



2 

1 c - 5 



G 



G 



07 



128 



j G\ 



05 



13 
128 



1G 

,G 2 



,G 



14 



16 



,G 



2:! 



25 
128' 



G 



34 



81 
128 



We can then compare the result for the boundary state 
found in this section with the specialization at c = 1/2 
of formula (28) by expressing the Virasoro modes in a 
standard way in terms of fermionic ones: 

Ln = \ X k : ^-k+nipk ■ , n G Z . (93) 
k e z+i 

We have verified the agreement for the first few descen- 
dants up to level 8 in the Virasoro modes of the two 
derivations for the Ising CFT. 

Further by comparing the amplitudes of boundary 
states we find: 



{B o 1^0-1/48^0) = 



-1/24 



13 



55 



1235 
2048' 



(94) 



10 



where \B^) as in cq. (87) and L is expressed in terms of 
fcrmions using eq. (93). Of course, we know on the other 
hand from the general case (53) that 



- 1/48 is;> 



-1/4 



(95) 



Matching (94) and (95) gives rise to intriguing combina- 
torial identities, which we have of course checked to high 
order, but are not able to prove. 



V. LATTICE MODELS 
A. Lattice discretization and numerics 

We propose now a lattice discretization of the bound- 
ary state and verify numerically the contribution of the 
first levels descendants of the identity in (28) by com- 
puting the scaling limit of scalar products in the lattice 
model. 

The lattice model we consider is the dense loop model 
based on the adjoint representation of the Tcmperley- 
Lieb algebra acting on TV (even) strands TLat(/3) (/? is 
the so-called weight of loops), and imposing free bound- 
ary conditions at both boundaries. At the critical point, 
the anisotropic version of the model is defined by the 
Hamiltonian: 



H 



JV-2 

E< 

i=0 



(96) 



where e; arc the TL generators. This Hamiltonian acts 
on link states |s)jvj reduced TL diagrams keeping track 
of connectivities of sites. The Hamiltonian can be put in 
a triangular block form, with j = 0, 1, . . . , TV, the num- 
ber of through lines, indexing each block [23]. On the 
lattice we define the transpose n(s\ of a state by turning 
it upside down and the bilinear form which we call loop 
scalar product n(si\s2)ni associating to the link states 
the diagram obtained by gluing |s2)iv with at(si|. The 
result is given by /3 n , with n the number of loops thus 
formed, and we get zero if we contract two strings. For 
instance we have 



(97) 



This is the usual bilinear form used for the TL algebra 
[24] . For generic j3 the loop scalar product is not positive 
definite. There is no reason to expect otherwise, as the 
theories we are dealing with are not unitary. The contin- 
uum limit of this loop model, when we parametrize the 
loop weight as (3 = 2 cos(7r/(p+l)) (p> 1 is a real param- 
eter), is a CFT with central charge c = 1 — 6/(p(p + 1)). 
Further the bilinear form introduced above flows in the 
continuum limit to the Virasoro bilinear form, as already 
observed in [16]. This will allow us to measure the bound- 
ary state on the lattice. We claim that the lattice state 



renormalizing in the scaling limit to our boundary state 
introduced above is the following link pattern 



\B°) 



N 



) 



TV- 



a\B°) 



(98) 



where in the continuum limit procedure we discard non- 
universal contributions, and introduce the proportional- 
ity constant a. Note that the lattice state is normalized 
according to the form jv(-|-)jv and that only the sector 
j = without through lines couples to this state. 

To check this relation we have computed numerically 
the scaling limits of lattice scalar products N(B°\k)N of 
the boundary state with the fe-th eigenvector in the sec- 
tor j = 0, ordered according to its energy (|0)jv is the 
ground state). The vectors |fc)_/v are normalized to 1 us- 
ing the loop scalar product. In the j = sector, the 
continuum theory has field content given generically by 
the quotient of Virasoro Verma modules V\,\jV\-\ (with 
the usual Kac table notation V r . s for Verma modules) 
[3], and the scaling limit \k) of |fc)jv will be given by 
a combination of Virasoro descendants at a level deter- 
mined by the energy of \k). We identify then obviously 
\k = 0) as the vacuum |0) and (due to normalization) 
|1) = y/2/cL- 2 \0), |2) = yT72cT:_3|0). Unfortunately 
we do not know a priori which combination of Virasoro 
descendants contributes to |fc), for k > 2. We however 
do know the number of descendants at each level, since 
the above quotient amounts to setting L_ 1 |0) = in the 
generic Verma module. 

Due to the anomaly at the corners, we expect the fol- 
lowing scaling form of the lattice scalar products: 



-log( N {B°\k) N ) 
with 



aoTV+ai log TV+a2 



f J_ 

TV TV 2 \N 3 
(99) 



(100) 



Note the difference between formula (99) and the one for 
the scaling in the case of periodic boundary states \B P ), 
for which there is no logarithmic term and excited states 
corresponding to descendant fields decouple. 

In the following we take anyway a\ as a free parameter 
and fit our numerical results with (99) where we drop 
the term 0(TV~ 3 ). (B°\k) can then be extracted from 
a-2 = — \og((B°\k) /a). As usual we cannot determine 
the sign of the scalar product, which however depends 
on an arbitrary phase of \k) N . The numerical results are 
presented in table I. 

The value of a\ and a are obtained from the scaling 
of (73° |0). We get good agreement of the numerical re- 
sult for ax with the CFT prediction when p = 3, and the 
discrepancy increases with p. The constant a clearly de- 
pends on p, but we have not been able to determine if it 
is a universal quantity or not. Simulations for the critical 
XY spin chain in a transverse field show that it is inde- 
pendent of the parameter r tuning from c = 1/2 to c = 1, 
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P 


numerics 


CFT 




3 


-0.06238 ± 0.00137 


-0.0625 


ai 


4 


-0.08659 ± 0.00214 


-0.0875 


5 


-0.09515 ± 0.00033 


-0.1 




OO 


-0.11706 ± 0.00762 


-0.125 




3 


0.97642 ± 0.00447 




Q 


4 


0.94890 ± 0.00679 




5 

oo 


0.94000 ± 0.00010 

n 00701 i f\ nnoo7 

0.88731 ± 0.00327 






3 


0.49994 ± 0.00229 


0.5 


(B«|l> 


4 


0.58652 ± 0.00422 


« 0.591608 


5 


0.61994 ± 0.00067 


« 0.632456 




oo 


0.65330 ± 0.01789 


» 0.707107 




3 







(B°|2) 


4 
5 

oo 











TABLE I. Numerical results for the loop model obtained by 
fitting the lattice scalar products with (99) for system sizes 
N = 8 -» 24. 



suggesting that it could be universal. Further it is shown 
in [14] that the ratio of two proportionality constants 
a for lattice boundary states corresponding to different 
boundary conditions is universal, and so we believe that 
a should contain a universal part. For computing (B°\k), 
we actually fit 



log 



N 



[B°\k) 



N 



N 



[B°\0) 



N 



cla ac. 
a2 + i + A^ + ^' (101) 



and drop first points. According to the identification 
|1) = y / 2/cL_2|0), the numbers {B°\l) in the column 
CFT are then yc/2. As for the comparison of ax, when 
p = 3 the agreement is very good, while it gets worse for 
greater values of p. |2) is a state with conformal dimen- 
sion 3, so that our formula of the boundary state predicts 
its decoupling. Remarkably, (B°\2) is exactly zero within 
numerical precision for sizes N > 10. For higher excited 
states numerical simulations do not give results accurate 
enough to allow to clearly identify the combinations of 
Virasoro descendants to which they correspond, so that 
we cannot compare further the prediction of our formula 
for the boundary state (28). Finally, we note that for 
p = 1,2 corresponding respectively to dense polymers 
(c = —2) and percolation (c = 0), the presence of null 
vectors in the module of the identity require smart tricks 
for computing scalar products, such as those employed 
in [16], and we have not studied numerically these cases. 



B. Ising chain 

In this section we give evidence that the fermionic co- 
herent state of eq. (87) is the universal part of the con- 
tinuum limit of a lattice state expressed in terms of Ising 
spins. Similar computations appear in [25]. D boundary 
conditions considered in the previous derivation corre- 
spond to fixed boundary conditions for the spins. As we 



have already observed, result (87) holds more generally 
for homogeneous boundary conditions, and here we will 
consider instead free boundary conditions (as done for 
the loop model in section VA). 

We consider the Hamiltonian limit of the 2D critical 
Ising model on a strip with free/free boundary conditions: 



/ N N-l \ 



(102) 



acting on an Hilbert space made of N (even) copies of 
the fundamental representation of su(2). This model 
has been extensively studied in the past (see for example 
[26]). After briefly reviewing the solution, we will com- 
pute the limit of lattice scalar products which we will 
compare with the CFT predictions. 

We introduce the fermions Cj,cj and perform the 
Jordan- Wigncr transformation (redefining of as —erf) 



of = I -24a, of = JJ(l-2ctc J -)(c i + cJ), (103) 
to obtain an expression of H quadratic in cj : 

N N-l 



C i C i+l "T" c i+l c i ! 

(104) 

which can readily be diagonalized through a canonical 
transformation 



N 

^±4 = E^( c * ±c 



Wc have 



/v 



with 

A fc = 2 sin 



in - 



fc=i 



2k - 1 



\2(2N + 1 
2 



/ 2k - 1 

COS 7T I I 



^/2N + 1 \2N+1 

(-r +1 — 



. (2k-\ 
sin iri 



V2N+1 \2N+1 



(105) 

(106) 

(107) 
(108) 
(109) 



Note that <fi are orthogonal matrices. The ground state 
of the Hamiltonian is 10)^, with ipk\0) N = for every k. 
When N is large, we have 



ea^^eMVl 

k=i fe=i ^ ' 



(110) 



If now we introduce the notation = ipj, for k > 0, 
and redefine tpk —> ^k-1/2, wc identify H with Lq, where 
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Lo is defined in terms of rcscaled fermions as in (93). 
We recover then that the system in the continuum limit 
is described by Ising CFT. For the free boundary condi- 
tions chosen the states propagating in the strip are the 
descendants of the identity (sector with an even num- 
ber of fermions) and the energy ip (sector with an odd 
number of fermions). 

Now recall that the Hamiltonian (102) is related to 
the 2D Ising model in the a x basis. Then to make the 
contact with the 2D model, we should rotate 90 degrees 
clockwise the spins of the chain in the x — z plane, that 
is, we define the free boundary state |SS)jv &s 



\B^)n = 



1= E 



| M ?.../4) = |t---t) . (Ill) 



-} 



We want to compute ^(B^\k) n, where |/c)at is the fc-th 
excited state. One way to deal with that is writing the 
projector onto |-B^)jv in terms of Cj's (the same trick was 
used for example in [27]) 



N 



\B^,)nn{B^\ - J 



i-oi 



N 

n 



eje,- 



(112) 



so that \N(B^\k)N\ 2 can be evaluated efficiently using 
Wick's theorem. Let us start dealing with the ground 
state. Rearranging the factors of cj and Cj in a convenient 
way, we have to compute 



N 



N 



N 



where (see also [28]) 



Gi 



Di 



-C,-. 



(44) 



G, 



C D 

-D -C 



Gji 



25, 



Gin — G^- 



and we have introduced the function 
G 



N 

fc=i 



(113) 

(114) 
(115) 

(116) 



Computing the Pfaffian from the square of the determi- 
nant and using properties of block matrices, we have 



N 



'M\k) 



AM 



det 

l<i,j<N 



Gi 



(117) 



The scalar product with the excited state ipkt ■ ■ ■ V'fes |0)jv 
can be computed by regarding this state as the vac- 
uum state of a new set of ip's, where ip^ and ipj, are 
interchanged for ki,...,k s [26]. For these A;'s we have 
Afe — ^ — Afc and 4>Zi — > — <fe, so that (117) holds also for 
ipki ' ' ' V'fcs |0)jy if we replace G^- with 



Gij [ki , . . . , k s ) 



^k^ij 

k unexc. 



+ E Wtr ( 118 ) 

k cxc. 



We have not been able to find a closed formula for 
the scalar products, but eq. (117) allows for very effi- 
cient numerical computations of the scalar products. We 
have fitted results for N = 2 — > 500 (dropping the first 
points for higher excited states), using the fit function 
(99). The value a\ = —1/16 is always found with a pre- 
cision of 10 -6 and for the proportionality constant a we 
find with a high precision a = 1.01937. For jv(_B^|fc) n we 
have the results in table II. We note that obviously the 





h k 


numerics 


CFT 


Wll) 


1/2 








{B$\2) 


3/2 










2 


0.499994 ± 0.000003 


1/2 




5/2 








W|5) 


3 








(w 


7/2 








(Bm 


4 


0.124995 ± 0.000003 


1/8 = 0.125 


{sm 


4 


0.624989 ± 0.000003 


5/8 = 0.625 




9/2 








(BZ\10) 


9/2 









TABLE II. Scalar product of \B1) and |fc) from finite size 
scaling (numerics) and comparison with CFT prediction, hk 
is the conformal dimension of the field | k) . 

sector of the energy -0 (hk half integer) decouples. The 
agreement is very good. It remains intriguing of course, 
that (95) does not seem to be easily obtainable from the 
explicit solution of the lattice model. This is in contrast 
with similar expressions for other geometries, such as the 
torus. 



VI. CONCLUSION 

This general exploration of the boundary states for the- 
ories defined on a segment will be used in our next paper 
to discuss in particular properties of geometrical prob- 
lems on a rectangle. Of course, the set-up can find many 
other applications as well, in condensed matter or string 
theory [9-11, 29]. Another field where our boundary state 
(28) might be relevant is the calculation of quantum in- 
formation quantities such as entanglement entropies or 
overlaps in Id or 2d systems. For example, in some 2d 
conformally invariant wave functions [27, 30], the calcula- 
tion of the entanglement entropy boils down to the study 
of expansions like (99). Some references [31, 32] have fo- 
cused on the universal logarithmic term a\ logL in (99), 
but we have shown that the next term ai contains also 
a universal piece, and it would be interesting to extend 
their analysis using our present results. 



A. Note added 

After the completion of this work, we became aware of 
a similar result published in [33-35], which is consistent 
with our main formula eq. (28). 
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